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ANALYSIS  OF  RESONANT  ORBITS 

by 
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Doreen  M.  C.  Walker 


,  SUMMARY 

. /*■  The  Earth's  gravitational  potential  is  usually  expressed  as  an  inf  in 

ite  series  of  tesseral  harmonics,  and  it  is  possible  to  evaluate  'lumped 
harmonics'  of  a  particular  order  m  by  analyses  of  resonant  satellite 
orbitsj-Jorbits  with  tracks  over  the  Earth  that  repeat  after  m  revolu¬ 
tions.  in  this  paper  we  review  results  on  30th-order  harmonics  from 
analyses  of  15th-order  resonance,  and  results  on  29th-  and  31st-order  har¬ 
monics  from  29:2  and  31:2  resonance. 

The  values  available  for  30th-order  lumped  harmonics  of  even  degree 
are  numerous  enough  to  allow  a  solution  for  individual  coefficients  of 
degree  up  to  40.  Th'  best-determined  coefficients  are  those  of  degree  30, 
namely 

1C,9e30,30  ’  -1*2  *  ‘-1  1(,9530,30  "  9’6  S  l'3 

(  The  standard  deviations  here  are  equivalent  to  1  cm  in  geoid  height. 

-For  the  29th-  and  31st-order  harmonics,  and  for  the  30th-order  har¬ 
monics  of  odd  degree,  there  are  not  enough  values  to  determine  individual 
coefficients,  but  the  lumped  values  from  particular  satellites  can  be  used 
for  'resonance  testing'  of  gravity  field  models,  particularly  the  Goddard 
Earth  Model  10B  (up  to  degree  36)  and  10C  (for  degree  greater  than  36). 

The  results  of  applying  these  tests  are  mixed. ^GEM  10B/C  emerges  well  for 
order  30,  with  sd  about  3  *  10"9;  for  order  31,7  the  GEM  10B  values  are 
probably  good  but  the  GEM  10C  values  are  probably  not;  for  order  29,  the 
test  is  indecisive. 
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1  INTRODUCTION 

A  satellite  orbit  experiences  6 : a  resonance  if  the  satellite  completes  g  revolu¬ 
tions  about  the  Earth  while  the  Earth  spins  a  times  relative  to  the  satellite’s  orbital 
plane.  The  satellite's  track  over  the  Earth  then  repeats  after  g  revolutions  in  approx¬ 
imately  a  days.  In  such  a  resonance  the  orbital  perturbations  caused  by  harmonics  of 
order  g  in  the  geopotential  have  the  same  effect  on  each  revolution,  so  that  they  build 
up  and  can  be  analysed  to  determine  accurate  values  of  'lumped  harmonics'  of  order  g  . 

The  lumped  harmonics  are  linear  sums  of  individual  harmonic  coefficients,  and  these 
individual  coefficients  can  therefore  be  determined  if  values  of  lumped  harmonics  can  be 
obtained  from  a  number  of  satellites  at  different  inclinations  to  the  equator. 

The  15:1  resonance,  or  15th-order  resonance,  when  the  track  repeats  after  15  revolu¬ 
tions  in  one  day,  has  proved  the  most  useful  in  practice,  because  the  height  (500-600  km) 
is  such  that  air  drag  brings  the  orbit  through  resonance  quite  slowly  -  but  not  so  slowly 
that  a  50-year  wait  is  needed.  An  orbit  that  stays  near  resonance  for  about  a  year  is 

ideal,  and  analysis  of  such  orbits  has  given  values  of  15th-order  lumped  harmonics  accurate 
1  2 

to  l at  best  *  .  Recently,  23  satellites  which  experienced  15th-order  resonance  were 
3 

used  to  determine  individual  harmonics  of  order  15  and  degree  15,  16,  17,  ...,  35.  Also 
values  of  individual  14th-order  coefficients  have  been  obtained  from  analysis  of  1 1 

4 

resonant  orbits  . 

The  2-day  resonances,  such  as  29:2  and  31:2,  are  much  weaker  than  the  14:1  or  15:1 

resonances,  but  the  effects  are  often  detectable.  The  first  such  analysis,  of  29:2 

resonance  by  Walker"*  in  1976,  has  been  followed  by  several  studies  of  29:2  and  31:2 
3-9 

resonance  ,  though  the  accuracies  achieved  have  not  been  better  than  about  15Z. 

Any  15th-order  resonant  orbit  which  stayB  near  resonance  for  a  year  or  more  is  also 
appreciably  affected  by  harmonics  of  order  30,  and  the  analysis  of  the  variations  yields 
values  of  30th~order  lumped  harmonics  as  a  by-product.  Again  their  accuracy  has  not 
usually  been  better  than  about  20X,  except  for  one  satellite,  where  71  was  achieved. 

The  aim  of  this  paper  is  to  gather  together  existing  results  on  lumped  harmonics  of 
order  29,  30  and  31  and  to  see  whether  it  ip  possible  to  produce  tentative  solutions  for 
individual  coefficients.  The  30th-order  harmonics  prove  to  be  the  most  satisfactory  and 
are  discussed  first  in  section  3,  after  a  brief  r£sum£  of  the  theory  in  section  2. 

Sections  4  and  5  are  devoted  to  the  29th  and  3!st-order  coefficients. 

2  THEORY 

The  longitude-dependent  part-  of  the  geopotential  at  an  exterior  point  (r,9,X)  can 
be  written  in  normalized  form*^  as 

r  Yj  Yj  ©  P”(C08  6){Sm  cos  mX  +  §*m  8in  mXH*  ’  U) 

fc-2  m-1 

where  r  is  the  distance  from  the  Earth's  centre,  6  is  co-latitude,  X  is  longitude 

.  3  2 

(positive  to  the  east),  u  is  the  gravitational  constant  for  the  Earth  (398600  km  / s  ) 


n|;iUJii  JU.iuom  ■".BJig>n^i»iy»-.HM.i'"!|.T;i»^!y»nyi<v|  .ji"  i’TOHJW|i|Jft!V  “i,i,,Tri  rv^-r- TT* * PiS^f T^nyfW! -  ’t*-w 
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and  R  ia  the  Garth's  equatorial  radius  (6378.1  km).  The  P®(cos  6)  are  the  associated 

Legendre  functions  of  order  m  and  degree  l  ,  and  C„  and  S„  are  the  normalized 

im  am  _ 

tesseral  harmonic  coefficients.  The  normalizing  factor  N^m  is  given  by 


N2  .  2(2a  4  D(a  -  m)l 
am  ( a  +  m) 1 


(2) 


A  satellite  experiences 
angle  defined  by 


6:a  resonance  when  i  •  0  ,  where 

$  -  o(o)  +M)  +  8  (ft  -  v) 


is  the  resonance 


Here  to  is  the  argument  of  perigee,  M  the  mean  anomaly,  (1  the  right  ascension  of  the 

» 

node,  and  v  the  sidereal  angle. 


Near  resonance  the  main  changes  in  inclination  1  due  to  the  relevant  C„  and 

_  am 

S.  take  the  form  of  terms  in  the  sine  and  cosine  of  4>,  2$ .  if  the  orbital  eccen- 

x.m  j  j 

tricity  e  is  small  .  The  multiplying  coefficients  of  these  terms,  which  we  seek  to 

_0,k 

evaluate  by  analysis  of  the  resonance,  are  the  lumped  harmonics,  written  as  C  and 
_0,k  m 

,  which  are  linear  sums  of  individual  coefficients  and  having  values 

of  a  that  increase  i  steps  of  2  with  a  %  m  and  (a  -k)  even.  The  index  k  is  equal 

to  yoi  where  y  »  1,  2,  3 .  and  the  terms  with  y  ~  1  are  usually  dominant. 


The  variation  in  eccentricity  e  is  controlled  mainly  by  terms  in  the  sine  and 

cosine  of  (*-q<o)  ,  the  terms  with  q  *  ±1  being  dominant.  Analysis  of  the  variations 
.  .  -q,k  -q»x 

in  eccentricity  yields  values  of  lumped  harmonics  Cm  and  ,  where  k  -  ya-q  : 

these  lumped  harmonics  are  also  linear  sums  of  individual  coefficients  C.  and  S„ 

am  am 

The  full  theory  has  been  given  previously  and  will  not  be  repeated  here.  Explicit 
forms  for  the  variations  in  i  and  e  at  29:2  and  31:2  resonance  can  be  found  in  Refs  12 
and  6.  Explicit  forms  for  the  30th-order  terms  at  15th-order  resonance  appear  in  Ref  2. 


3  HARMONICS  OF  ORDER  30 


3. 1  Form  of  the  even-degree  lumped  harmonics 


The  variation  of  inclination  at  15th-order  resonance  depends  primarily  on  the  15  th- 

order  harmonics  (the  terms  with  y  -  1  as  defined  in  section  2),  but  there  is  also  a 

.  .  2 
contribution  from  30th-order  terms  (y  ■  2) ,  and  the  main  30th-order  term  in  di/dt  is 


2n 


sin 


/r\30  _  (-0,2 

I  (a/  (15  “  c08  l)F30, 30,  Ur"  sin  2* 


_0,2 


30 


15 


-  S 


30 


cos  2$ 


15 


(3) 


i 

j 

A 


<3 

1 


Sr. 


where  n  (**  M)  is  the  mean  motion,  a  is  the  semi  major  axis,  —  to+M  +  15(0  —  v)  is 

the  15th-order  resonance  angle,  and  F  n  ..  is  Allan's  normalized  inclination 

j  j  >  1  4 

function  .  The  lumped  coefficients  are  given  by 


C  +  o0’2  C  +  Q^»2  r  + 

30,30  ”32  j2,30  v34  ^34,30 


(4) 


ja^YiiTfiirn' 


with  the  sane  equation  for  S  on  replacing  C  by  S  throughout.  The  Q-factors  are 
products  of  inclination  functions  F  and  powers  of  (R/a) ,  the  rth  Q-factor  in  equation 
(4)  being  given  by 


"‘30+2r 


(-# 


30+2r,30> 14+1/30,30,14 


In  each  analysis  of  a  15th-order  resonant  satellite  there  is  a  possibility  of 
extracting  determinate  values  for  these  even-degree  lumped  30th-order  coefficients,  though 
in  practice  only  the  best  resonances  yield  such  values.  If  values  of  lumped  coefficients 
can  be  obtained  at  enough  different  inclinations,  it  should  be  possible  to  solve  equations 


(4)  for  the  individual  harmonics  C^30  «nu  *£>30  v- 
3.2  Observational  values  of  even-degree  lumped  harmonica 

In  practice  we  have  found  only  seven  analyses  of  15th-order  resonance  giving  values 
of  30th-order  coefficients  that  we  regard  as  reliable.  They  are  listed  in  Table  l  and 
briefly  described  below. 

The  results  for  1974-34A,  Intercosmos  11,  are  from  the  analysis  of  129  weekly  US 
Navy  orbits  at  dates  between  June  1975  and  November  1977  (Ref  2).  The  analysis  yielded 
excellent  values  of  15th-order  harmonics,  but  the  two  30th-order  coefficients  are  small 
and  scarcely  greater  than  their  standard  deviation.  As  the  expression  (3)  shows,  however, 
the  effect  of  the  harmonics  is  best  indicated  by  multiplying  by  *"30,30, 14  *  w^ch  vari-es 
greatly  with  inclination.  Table  1  also  shows  the  values  after  such  multiplication  and, 
judged  in  this  light,  the  sd  for  1974-34A  emerges  as  the  second  best  among  the  seven 
satellites;  the  values  are  of  the  same  order  as  the  sd  only  because  the  values  happen  to 

be  small. 

The  results  for  1963-24B,  Tiros  7  rocket,  are  from  analysis  of  112  weekly  US  Navy 
orbits  ac  dates  between  November  1975  and  December  1977  (RefQ2).  Again  the  values  of  the 
15th-order  coefficients  were  excellent,  but  the  value  for  C3’  happens  to  be  small,  as 
is  seen  from  the  entry  in  Table  1  for  ^30,30, 14C30  ’ 

The  third  and  fourth  satellites,  1970-1UA  and  1971-13B,  were  in j almost  identical 
orbits  but  exhibited  variations  in  inclination  of  quite  different  form  ’  because 


( l  even) . 


resonance  occurred  at  different  values  cf  *,5  .  The  excellent  agreement  between  the 
values  of  C®’2  and  S°J2  obtained  for  the  two  satellites  gives  confidence  in  their 

reliability. 

The  fifth  satellite  is  1967-42At  Ariel  3,  which  has  served  as  the  standard  satellite 
in  Gooding's  development  of  the  computer  programs  THROE  and  SIMRES  used  m  resonance 
analysis.  The  values  in  Table  1  are  from  a  SIMRES  fitting  of  inclination  and  eccentricity 
on  281  orbits  at  dates  between  May  1967  and  August  1969,  with  ten  coefficients  . 

The  sixth  satellite  is  1971-54A,  for  which  269  weekly  US  Navy  orbits  over  five  years 
were  analysed1.  The  results  were  excellent  and  give  the  most  accurate  values  of  lumped 
30th-order  coefficients  yet  obtained. 


and  S°’2  obtained  for  the  two  satellites  gives  confidence  in  their 

in 
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The  final  orbit,  that  of  1964-S2B,  Nimbus  1  rockat,  was  datsrminad  by  Hiller'6  at 
25  epochs  between  March  and  September  1970,  from  2000  observations  including  many  Hewitt 
camera  plates.  Although  the  individual  orbits  are  more  accurate  than  those  of  the  other 
six  satellites,  only  41  orbits  were  used  in  the  resonance  analysis,  so  the  reliability 
of  the  values  of  the  lumped  coefficients  cannot  be  guaranteed. 


Table  1 

.  .  0  2 

Values  of  even-degree  lumped  harmonics  (C,S)^  for  the  seven  satellites 


No. 

Satellite 

i 

(deg) 

10  C30 

O-0.2 

10  S30 

*30,30,14 

o-  -0.2 

10  *30,30, I4C30 

9-  -0,2 
10  F30,30,I4S30 

1 

I974-34A 

50.64 

597  ±  558 

679  ±  651 

0.000952 

0.57  ±  0.53 

0.65  1  0.62 

2 

I963-24B 

58.20 

46  i  106 

253  1  88 

0.01176 

0.54  1  1.25 

2.98  i  1.03 

3 

1970-11 IA 

74.00 

19.2  ±  4.9 

4.1  t  4.4 

0.2579 

4.95  1  1.26 

1.06  1  1.13 

4 

197I-I3B 

74.05 

27.1  ±  5.5 

6.0  1  3.3 

0.2594 

7.03  ±  1.43 

1.56  1  0.86 

5 

I967-42A 

80.17 

-9.1  ±  4.6 

-5.0  t  5.5* 

0.4340 

-3.95  t  2.00 

-2.17  ±  2.39* 

6 

197I-54A 

90.21 

-8.2  t  0.6 

11. 1  10.8 

0.4755 

-3.90  ±  0.29 

5.28  i  0.38 

7 

I964-52B 

98.68 

22.8  i  7.9 

38.0  i  10.3* 

0.2502 

5.70  1  1.98 

9.51  i  2.58* 

*  These  standard  deviations  were  doubled  in  the  solutions 


3.3  The  equations  to  be  solved 

We  have  seven  equations  of  the  form  (4)  for  C  ,  and  seven  for  S  .  The  values  of 
the  Q  coefficients  for  the  seven  satellites,  evaluated  with  Gooding’s  computer  program 
PROF,  are  given  in  Table  2  on  page  7. 


In  addition,  following  the  procedure  that  has  proved  successful  in  the  past,  we  add 
constraint  equations  of  the  form 


E».»  *  0  4  l0'5/l2 
§M0  ‘  0  *  l0'5/*2 


(6) 


where  l  -  30,  32,  34,  ...  These  equations  express  the  expectation10  that  the  order  of 

.  .  .  ,  .  -S  2 

magnitude  of  the  individual  coefficients  of  degree  l  is  10  /t,  for  30  «  i  4  50  , 

as  is  confirmed  in  a  general  way  by  the  Goddard  Earth  Model  10C  (Ref  17), 

Thus  we  solve  7  +  N  equations  by  least  squares  for  N  harmonics,  the  optimum 

value  of  N  being  selected  empirically.  As  in  our  recent  solution  for  15th-order 
.  3 

harmonics  ,  we  are  prepared  to  consider  a  relaxation  of  some  standard  deviations,  if 
necessary. 


3.4  The  solutions  for  individual  coefficients  of  even  degree 


When  the  equations  were  solved,  there  were  no  problems  with  the  C-values,  but  two 
of  the  S -values,  for  1967-42A  and  1964-52B,  had  large  weignted  residuals,  for  all  values 
of  N  up  to  6.  The  standard  deviations  for  these  two  satellites  were  therefore  doubled 
as  indicated  in  Table  1 . 


I 
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Table  2 

0  2  0  2  0  2 

Value*  of  Q^’  »  Q-'i,  »  »•••  for  the  seven  satellites 


Satellite 

Q0*2 

^32 

Q0*2 

w34 

Q°*2 

y36 

Q0'2 

y38 

p0,2 

w40 

Q0*2 

^42 

1974-34A 

-11.9 

52.9 

-131.8 

200.3 

-171.5 

34.5 

1963-24B 

-7.61 

19.01 

-20.11 

2.71 

11.18 

-2.92 

1970-1  1A 

-1.121 

-0.806 

-0.128 

0.319 

0.433 

0.314 

1971—1 3B 

-1.107 

-0.807 

-0.140 

0.308 

0.429 

0.317 

1967-42A 

0.155 

-0.161 

-0.281 

-0.295 

-0.252 

-a- 

CO 

? 

1971-54A 

0.430 

0.213 

0.100 

0.036 

-0.000 

-0.020 

1964-52B 

-1.113 

-0.847 

-G.231 

0.127 

0.313 

0.333 

When  the  7  +  N  equations  were  solved  for  N  coefficients,  for  2  <  N  <  6  ,  the 
values  of  the  measure  of  fit  e  were  as  follows: 

N  2  3  4  5  6 

C  equations  1.08  1.05  1.04  0.92  0.89 

S  equations  1.50  0.96  0.95  0.94  0.92 

o 

As  usual,  e  is  the  sum  of  squares  of  the  weighted  residuals,  divided  by  the  number  of 
degrees  of  freedom;  and  the  weighted  residual  is  the  residual  for  each  lumped  coefficient 
divided  by  the  standard  deviation  for  that  coefficient,  as  given  in  Table  1. 

Thus  it  appears  that  any  value  of  N  between  3  and  6  wil-  give  a  'satisfactory' 
solution  with  c  near  1.  However,  Table  2  shows  that  1974-34A,  the  second  most  accurate 
of  the  satellites,  has  its  largest  Q  coefficients  for  i  -  38  and  t  ■  40  .  So  the 
6-coefficient  solution  (Jl  ■  30,  32,  ....  40)  is  required  to  ensure  that  1974-34A  plays  its 
proper  role. 

Table  3  gives  the  6-  and  4-coefficient  solutions,  together  with  the  values  from  the 
Goddard  Earth  Model  10B  (Ref  18),  which  only  goes  to  degree  36  and  so  should  be  compared 
with  the  4-coefficient  solution,  and  also  values  for  l  ■  38  and  l  »  40  from  GEM  10C 
(in  brackets) 

Table  3 

Values  of  8n<*  S£  30  8iv"n  by  our  6-  and  4-coefficient 

solutions  and  by  GEM  10B  (and  GEM  10CT 


1 

,0\xi 

6-coefficient 

4-coefficient 

GEM  10B 
(GEM  10C) 

6-coefficient 

4-coefficient 

GEM  I0B 
(GEM  IOC) 

30 

-1.2  ±  1.1 

-0.4  ±  1.0 

-5.2 

9.6  ±  1.3 

10.1  ±  0.9 

11.1 

32 

-14.5  ±  4.3 

-16.2  t  4.1 

-0.6 

0.1  +  4.4 

-1.4  ±  3.6 

-0.2 

34 

-4.7  ±  4.4 

-3.6  ±  4.8 

-11.9 

7.2  ±  4.6 

7.9  +  4.4 

1.2 

36 

0.7  ±  4.4 

-1.4  ±  4.2 

-3.9 

-2.4  ±  5.1 

-2.3  ±  4.3 

-0.9 

38 

5.7  ±  3.6 

(0.6) 

2.0  ±  4.5 

(2.6) 

40 

2.5  ±  3.7 

(1.5) 

2.4  i  4.1 

(-3.0) 

Table  3  shows  that  chs  solutions  srs  stable!  going  from  the  4-  to  the  6-cotff icient 
solution  does  not  change  any  value  by  as  much  as  its  sd.  This  is  largely  because 
I974-34A  dominates  for  i  ■  38  and  40,  and,  since  it  has  small  lumped  coefficients, 
I974-34A  tends  to  produce  small  38th-  and  40th-degree  coefficients. 

3 

The  15th-order  coefficients  in  GEM  I  OB  were  judged  to  be  probably  accurate  to  about 
3  x  10  If  the  same  figure  is  used  for  the  30th-order  coefficients  of  GEM  I0B,  Table  3 
shows  that  the  mean  difference  between  a  GEM  10B  value  and  the  corresponding  value  in  our 
4-coefficient  solution  is  0.8  * (the  sum  of  the  sd) .  Since  the  GEM  10B  values  and  ours 
are  believed  to  be  completely  independent,  this  agreement  suggests  that  a  standard  devia- 

.Q 

tion  near  3  *  10  may  be  appropriate  for  the  30th-order  GEM  10B  coefficients.  The  four 

-9 

GEM  IOC  values  given  in  Table  3  are  also  consistent  with  an  assumed  accuracy  of  3  «  10 
though  a  comparison  of  only  4  out  of  32000  coefficients  can  scarcely  be  regarded  as  a 
useful  test. 


However,  the  main  conclusion  from  Table  3  must  be  that  the  values  in  our  ac.'dtion 
are  poorly  determined,  except  for  *  ■  30  .  This  poor  accuracy  is  to  be  expected  when 
using  only  seven  satellites,  of  which  three  (1970-1 1 1A,  1971—1 3B  and  1964-52B)  have  very 
similar  values  for  their  Q-factors  (see  Table  2). 

The  weighted  residuals  in  the  6-coefficient  solution,  for  the  seven  satellite 
equations  and  the  six  constraint  equations,  are  givui  in  Table  4.  It  will  be  seen  that 
the  largest  weighted  residual  is  1.60. 

Table  4 

Weighted  residuals  in  the  13  equations  used  in  the  6-coefficient  solution 


Satellite  equations 

Constraint  equations 

_0, 2 

_0, 2 

Satellite 

C30 

S30 

Degree  * 

S,30 

**,30 

1974-34A 

0. 10 

0,01 

30 

0.1 1 

-0.87 

1963-24B 

-0.02 

0.30 

32 

1 .48 

-0.01 

1970-1 1 IA 

-0.50 

-0.36 

34 

0.54 

-0.83 

1971—1 3B 

1 .04 

0.09 

36 

-0.10 

0.31 

1967-4 2A 

-0.84 

-1.18 

33 

-0.83 

-0.28 

1971-54A 

0.03 

0.13 

40 

-0.40 

-0.40 

1964-52B 

0.33 

1.60 

3.5  The  variation  of  even-degree  lumped  harmonics  with  inclination 

Our  basic  data  are  the  seven  values  of  lumped  harmonics,  and  the  most  useful  compari 

sons  can  be  made  by  considering  their  variation  with  inclination  after  multiplication  by 

F30  30  14  *  1  8hows  the  values  of  the  lumped  harmonics  with  their  sd,  together  with 

the  variations  given  by  the  6-coefficient  solution  (unbroken  line),  the  4-coefficient 

0  2 

(dot-dash  line)  and  GEM  10B  (broken  line).  The  values  of  S^q  for  1967-42A  and  1964-52B 
are  shown  in  Fig  1  with  their  standard  deviations  doubled,  and  it  is  clear  that  the 
relaxation  was  justified. 


It  is  worth  noting  that  the  curves  of  Fig  1  are  nearly  symmetrical  about  i  ■  86°  , 
so  that  the  98°  satellite  (I964-52B)  competes  with  the  two  at  74°,  while  the  80°  satellite 
(1967-42A)  competes  with  that  at  90°. 

Fig  l  shows  that  our  6-coefficient  solution  is  only  slightly  better  than  thi  4- 


coefficient  solution,  and,  if  this  were  the  only  criterion,  the  latter  would  be 
recoomended .  However,  as  already  explained,  six  coefficients  are  required  to  take  account 
of  1974-34A,  and,  since  they  do  not  cause  instability,  the  6-coefficient  solution  is  to  be 
preferred. 

Another  notable  feature  of  Fig  1  is  the  good  performance  of  GEM  10B:  although  there 
are  some  discrepancies,  GEM  I  OB  generally  gives  a  variation  that  can  fairly  be  called 
quite  realistic. 


Review  of  values  of 


'30,30 


30,30 


Comparisons  with  previous  evaluations  of  individual  coefficients  are  not  likely  to 
be  very  illuminating,  because  of  the  general  inaccuracy.  But  comparisons  are  worthwhile 
for  C.jq  20  and  S3Q  30  *  which  are  well  determined  here. 

Table  5  lists  values  of  these  coefficients  from  various  sources.  The  first,  from  the 

European  Gravity  Field  Model  GRIM  2  (Ref  19),  would  not  be  expected  to  be  accurate,  because 

it  is  the  final  coefficient  in  a  30  *  30  array  and  is  likely  to  carry  a  burden  of  error 

from  the  undetermined  coefficients  of  higher  order  and  degree.  The  second  entry  is  from 

Rapp's  analysis  of  5°  *  5°  terrestrial  gravi  ,  measurements20,  which  goes  to  order  and 

degree  52  and  should  therefore  avoid  this  source  of  error.  The  same  applies  to  the  next 

entry,  GEM  10B,  which  is  complete  to  order  and  degree  36.  (GEM  10C  is  not  listed  because 

22 

it  is  the  same  as  GEM  10B  to  order  and  degree  36.  The  gravity  field  derived  by 
Gaposchkin  (1980)  goes  to  degree  and  order  30,  but  has  no  (30,30)  coefficient.)  All  these 
values  are  independent  of  ours . 

21 

The  values  obtained  by  Kosteleck^  and  Klokodnfk  ,  however,  utilize  some  of  our 
resonances  and  are  therefore  not  independent  of  our  values.  Kosteleck#  and  KlokoSnik 
took  pairs  of  lumped  coefficients  from  ten  satellites,  of  which  eight  pairs  were  obtained 
by  us.  Of  these  eight,  there  are  three  for  which  we  no  longer  regard  the  30th-order 
harmonics  as  being  reliably  separated  from  the  15th-order,  as  explained  in  Ref  3.  These 
three  are  Saturn  SA5,  1964-05A  (i  -  31°),  which  was  omitted  even  from  our  15th-order 
analysis;  Ariel  1,  1962—1 5A  (i  ■  53.8°),  and  Cosmos  72,  1965-53B  (i  ■  56.0°).  Of  the  other 
five  satellites,  we  use  three  here  with  the  same  values  as  Kosteleck^  and  Klokocnik;  and 
we  use  two,  1967-42A  and  1971-54A,  with  new  values.  Kosteleck$r  and  KlokoSnik  also  took 
two  values  from  their  own  analyses.  We  have  not  used  these  because  they  are  both  at 
inclinations  where  we  have  more  accurate  values  (50.6  and  74.0  ).  To  summarize, 

Kosteleckjr  and  KlokoJnik  use  the  same  values  as  we  do  for  three  satellites,  together  with 
three  of  our  values  we  have  now  discarded,  two  of  our  values  we  have  now  revised,  and  two 
of  their  own. 
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Table  3 

Values  of  C303Q  and  S3Q>3() 


Date 

Source,  and  type  of  data 

10  So, 30 

10  ^30,30 

1976 

GRIM  2  (Ref  19),  comprehensive 

14.2 

5.8 

1977 

Rapp  (Ref  20) ,  terrestrial  gravity 

1.4  ±  5.2 

12.6  t  4.9 

1978 

GEM  10B  (Ref  18),  comprehensive 

-5.2  ±  3? 

11.1  ±  3? 

1°79 

/Kosteleck?  i  (Ref  21  1  /solution  5 

(and  KlokoCnik f  (resonance/  (solution  7 

-6.0  ±  0.9 
-8.2  ±  5.7 

11.0  ±  1.8 
8.4  ±  5.8 

1981 

-  „  ,  .  /4-coefficient 

Our  solution,  resonance,  (6_coe£ficient 

-0.4  t  1.0 
-1.2  ±  1.1 

10.1  t  0.9 
9.6  ±  1.3 

How  do  the  values  in  Table  5  stand  up  to  comparison?  The  answer  is  "surprisingly 
well",  when  the  earliest  (GRIM  2)  is  ignored.  If  the  values  from  our  solutions  are 
accepted  as  correct,  averaged,  and  given  a  'safety  factor*  of  2  on  the  ad,  we  can  say  that 
the  C  and  S  values  are  -112  and  10  ±  2  respectively.  These  are  consistent  with 
nearly  all  the  previous  values  and  there  ia  remarkable  unanimity  on  the  value  of 

9_ 

10  S3q  3q  ,  all  six  values  being  between  8.4  and  12.6. 

3.7  Another  look  at  the  two  rejected  satellites 

We  rejected  Ariel  1  and  Cosmos  72  in  our  analysis,  because  we  preferred  a  fitting 
of  the  15th-order  resonance  which  had  no  30th-order  terms.  In  other  words,  we  judged 
that  the  30th-order  terms  did  not  contribute  usefully  to  the  fitting.  However,  this  does 
not  necessarily  mean  that  they  are  wrong,  and  it  is  of  interest  to  see  how  they  fit. 

Ariel  1  gives^  lO^FC  -  -4.0  ±  3.4  and  lO^FS  ■  2.2  ±  3.8  at  i  »  53.8°  ,  where 

the  values  from  the  unbroken  curve  of  Fig  1  are  0.8  and  1.6  respectively.  The  'residuals' 

24  9 _ 

are  therefore  -1 .4o  and  0.2o  respectively.  Cosmos  72  gives  10  FC  »  -4,1  ±  2.1  and 

9 _ 

10  FS  »  -1.7  ±  2.2  at  i  ■  56.0  ,  where  the  values  from  the  unbroken  curve  of  Fig  1  are 

0.7  and  2.2  respectively.  The  'residuals'  are  therefore  -2.3c  and  -1.8o  respectively. 

Thus  our  decision  to  exclude  these  two  satellites  seems  reasonable,  even  with  the 
benefit  of  hindsight.  It  would  be  possible  to  strengthen  the  S-solutions  by  including 
the  value  from  Ariel  1;  but  this  would  be  capricious,  since  we  had  previously  rejected 
the  solution  in  which  it  occurs. 

3.8  Lumped  harmonics  of  odd  degree 

So  far  we  have  only  discussed  the  coefficients  of  even  degree  (30,  32,  34,  ...), 
obtained  from  analysis  of  inclination.  Lumped  coefficients  of  odd  degree  (31,  33,  j,,  ...) 
can  in  principle  be  obtained  from  analysis  of  eccentricity  at  15th-order  reso.iance 
(explicit  forms  are  given  in  Table  4  of  Ref  2) .  But  in  practice  the  effects  of  these 
30th-order  harmonics  are  small,  and  only  one  satellite,  1971-54A,  has  yielded  determinate 
values  of  lumped  coefficients  of  odd  degree*. 
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Fig  2  shows  ths  vsluss  of  ths  lumpsd  coefficients  (C,§)  ;  and  (C,S)“q*3  detev- 
minsd  from  197I-54A  and  also  the  variations  of  thasa  coafficiants  with  inclination  piven 
by  GEM  10B  and  GEM  IOC .  Sines  coafficiants  of  degree  higher  than  36  contribute  signif¬ 
icantly  to  the  lumped  coefficients  at  most  inclinations!  and  GEM  I0B  does  not  go  beyond 
dagraa  36 ,  there  are  considerable  differences  between  the  GEM  I0B  and  GEM  IOC  curves. 

Both  curves  agree  fairly  well  with  the  S-coefficients  from  I9’l-54A,  but  the  C-coafficients 
of  GEM  appear  to  be  too  small. 

4  HARMONICS  OF  ORDER  29 


4.1  Form  of  the  even-degree  lumped  harmonics 

At  the  29:2  resonance  the  variation  of  i  it  given  by 
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Ai  n  /n\30  -  l-0*2  -0*2  I 

dt  *  sinT  (a)  (29  "  2  c0#  i)FjO,29,14jS29  8ln  <J,29  +  C29  co*  *29f 


^  terms  in  e ,  e  ,  . . . , 


(7) 


where  $2^  *  2(u+M)  +  29(0  -  v)  (8) 

and  the  lumped  coefficients  are  given,  by 


'30,29 


F32,29,15 
F30 ,29,14 


F30,29,14 


•  •  • 


(9) 


with  the  same  equation  for  S  on  replacing  C  by  S  throughout.  Equation  (9)  shows 
that  the  variation  of  inclination  depends  on  the  harmonics  of  even  degree. 

4.2  Observational  values  of  even-degree  lumped  harmonics 

There  are  four  satellites  for  which  the  orbital  inclination  has  been  analysed 
during  passage  through  29:2  resonance.  These  are  discussed  in  turn  below,  in  order  of 
increasing  inclination. 

The  phenomenon  of  29:2  resonance  was  first  recognised  and  analysed  for  Ariel  1, 
1962-15A,  by  Walker3.  The  analysis  showed  the  reality  of  the  phenomenon,  and  the  possi¬ 
bility  of  evaluating  lumped  harmonics;  but  the  data  were  not  of  high  accuracy  and  the 
values  of  lumped  coefficients  obtained,  given  in  Table  6,  were  inevitably  rather 
inaccurate. 

9 

The  second  satellite  is  1976-62E,  Cosmos  837  rocket,  recently  analysed  by  Hiller  . 

This  is  the  best  example  of  29:2  resonance  so  far  studied.  But  it  is  still  not  ideal 

because  the  eccentricity  was  appreciable  (0.04)  and  the  e  terms  in  equation  (7)  should 

strictly  have  been  included.  This  was  not  practicable  because  the  inclination  was  close 

to  63.4°  where  <J  ■  0  ,  and  consequently  the  e  terms,  which  depend  on  t  u)  , 

•  cos  "  n  a 

were  strongly  correlated  with  the  main  terms,  in  ,  *,Q  . 

8  in  *7 
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The  third  and  fourth  satellites  in  Table  6,  1971-106A  and  1971-18B ,  were  both  of 

high  drag.  The  changes  in  inclination  at  resonance  were  therefore  very  small  and  the 

1 2  23 

values  obtained  for  the  lumped  coefficients  have  large  standard  deviations  ’ 

If  the  early  terms  were  dominant  in  equation  (9),  it  might  be  possible  to  use  the 
values  of  the  lumped  coefficients  from  these  four  satellites  to  solve  for  two  or  three 
individual  coefficients;  but  at  inclinations  between  50°  and  65°  the  higher-degree  terms 
are  important,  and  often  dominate,  so  a  solution  is  not  practicable.  The  best  that  can 
be  done  is  to  compare  the  values  with  those  from  GEM  IOC. 

Table  6 

-  -  0  2 

Values  of  lumped  harmonics  (C,S)29  from  the  four  satellites 


Satellite 

i 

(deg) 

q  0,2 

io  c29 

q.0,2 

i°  s29 

F30,29, 14 

o-  -0.2 

10  F30,29, 14C29 

9-  -0,2 
10  F30,29,14S29 

1962-15A 

53.8 

-1000  ±  500 

1100  ±  500 

0.01526 

-15.3  ±  7.6 

16.8  ±  7.6 

1976-62E 

62.7 

-10  ±  15 

-76  ±  12 

0.1286 

-1.3  ±  1 .9 

-9.8  ±  1.5 

1971-106A 

65.7 

-90  ±  74 

-127  ±  39 

0.2089 

-18.8  ±  15.5 

-26.5  +  8.1 

197 1  —  1  SB 

69.9 

50  ±  20 

-160  ±  70 

0.3372 

16.9  ±  6.7 

-54.0  ±  23.6 

4.3  The  variation  of  even-degree  lumped  harmonics  with  inclination 

.  -  -  -  0  2 
Fig  3  shows  the  four  pairs  of  values  of  F3Q  29  ^(C.S)^  from  Table  6,  and  the 

curve  of  variation  with  inclination  as  given  by  GEM  10C,  with  coefficients  up  to  degree  44 
included.  The  comparison  is  inconclusive.  It  could  be  argued  that  most  of  the  observed 
values  are  numerically  too  large.  On  the  other  hand,  it  may  be  that  the  coefficients  of 
order  29  in  GEM  IOC  are  much  less  accurate  than  those  of  order  30:  certainly  the  S- 
coefficients  in  GEM  IOC  appear  to  be  too  small,  and  GEM  IOC  fails  to  fit  the  well  deter¬ 
mined  FS  value  of  -9.8  ±1.5  at  inclination  62.7°.  Also  Rapp’s  S-values20  are  much 
larger  and  give  a  variation  quite  different  from  the  GEM  10C  curve. 

Further  observational  result  are  needed  before  any  conclusions  can  be  drawn. 


Lumped  harmonics  of  odd  degree 


Lumped  coefficients  of  odd  degree  can,  in  principle,  be  determined  from  the  varia¬ 
tions  in  eccentricity  at  29:2  resonance,  but  in  practice  the  effects  are  small,  and  the 
values  so  far  obtained,  from  1962-1 5A  and  1976-62E,  cannot  be  regarded  as  well-established. 
It  is  not  worth  attempting  any  comparisons  until  better  observational  values  are  available: 
it  is  hoped  that  analysis  of  the  orbit  of  1968-40B,  which  is  currently  passing  through 
29:2  resonance,  will  provide  good  values. 
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5  HARMONICS  OF  ORDER  31 

5.1  Form  of  the  lumped  harmonics 

At  the  31:2  resonance  the  variation  of  i  is  given  by7 

..  ,R.32  _  (0,2  _0,2 

dt  "  linT  (a)  (31  "  2  cos  x^32,31 , 15/31  sin  *3!  +  Si  cos  *3 


where 


*3)  -  2 (w  +  M)  +  31(0  -v) 


+  terms  in  e,  e  ,  (10) 

(11) 


and  the  lumped  coefficients  are  given  by 
-®*2  -  Sa  Ti  l*  /p\2 


c32  -  lablLJfi  gtg  +  l^JLLdl  5 

32,31  ~  \fl/  34,31  n  \ fl  /  36, % 

32,31,15  32,31,15 


with  the  same  equation  for  S  on  replacing  C  by  S  throughout.  Equation  (12)  shows 
that  changes  in  inclination  again  yield  values  of  harmonics  of  even  degree. 

However  it  has  also  been  possible  to  analyse  the  variation  of  eccentricity  at  31:2 
resonance  for  three  satellites,  to  obtain  values  of  lumped  harmonics  of  odd  degree.  The 
variation  of  e  is  given  by7 


+  13F, 


,1,1 

_1>1 

'31  Sin<*31 

"  ->  -  S3. 

cos($3] 

-1,3 

-1 

,3 

'31  8in(*31 

+  “  531 

cos(4> 

+  terms  in  e,  e  ,  . . . , 


The  lumped  harmonic  coefficients  here  are  given  by 
>»•  _  >8F„  „  /n\2  _  19F, 


_-l  *3 

C„. 


l8^33,31 

ill  ( 

if  c  + 

19F 

35,31,17 

17*31,31 

\ 

,15 

,a/  33,31 

1  7F 

1 31 ,31 ,15 

14F 

'33,31 

ill  ( 

if  c  + 

1 5F 

3  35,31,16 

1 3F 

* J  31 ,31 

\ 

,14 

,a/  °33, 31 

1 3F 

*"31,31,14 

L  (!)  Ss, 


(!)  Sv 


and  similarly  for  S  on  replacing  C  by  S  throughout. 

5.2  Observational  values 

So  far,  four  orbits  have  been  analysed  at  31:2  resonance.  They  are  discussed 
overleaf  in  order  of  increasing  inclination. 
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The  variations  in  inclination  and  eccentricity  for  Skylab  I  rocket,  1973-27B,  at 
31:2  resonance  were  analysed  by  King-Hele^,  and  quite  well-determined  values  were  obtained 
for  all  six  of  the  lumped  coefficients.  These  values  are  given  in  Tables  7  and  8. 

The  orbit  of  Proton  4,  1968-103A,  was  analysed  by  Hiller  and  King-Hele^.  Since  the 
orbit  was  appreciably  eccentric,  it  was  found  necessary  to  include  the  e  terms  in 
equation  (12),  and  the  ’main’  coefficients  and  )  turned  out  to  be  very  small 

(they  are  given  in  square  brackets  in  Table  7).  A  simultaneous  analysis  of  inclination 
and  eccentricity  gave  the  values  for  the  other  four  coefficients  listed  in  Table  8. 

The  third  satellite,  China  2  rocket,  1971-18B,  suffered  high  drag  and  the  values 

.  .  25 

obtained  from  analysis  of  inclination  are  barely  determinate  .  They  are  included  in 
Table  7. 

8 

The  fourth  and  best  analysis  is  of  Samos  2,  1961  a 1,  by  Walker  .  A  simultaneous 
fitting  of  inclination  and  eccentricity  was  successful  in  giving  well-defined  values  of 
all  six  lumped  coefficients,  as  recorded  in  Tables  7  and  8. 

It  is  not  possible  to  solve  for  individual  coefficients,  because  the  numerical 
factors  multiplying  the  high-degree  terms  are  extremely  large  for  I973-27B  and  1968-103A. 
For  example^,  with  1973-27B,  has  a  maximum  numerical  factor  of  979.6  for  C^j  3)  » 

and  has  a  factor  363.1,  so  that  coefficients  would  need  to  be  evaluated  up  to 

degree  57,  and  this  is  impossible  with  only  six  equations. 

5.3 

Fig  4  shows  the  four  pairs  of  values  of  F32  3]  j 5 (C ,S) 3 *  from  Table  7  and  also 

the  variation  of  these  lumped  harmonics  with  inclination  as  given  by  GEM  IOC  up  to 
degree  46. 

•  ——II  ——-|Q 

Fig  5  shows  the  values  of  the  four  lumped  harmonics  (C ,S) and  (C,S)3J'  , 

multiplied  by  the  appropriate  F  ,  for  each  of  the  three  satellites  of  Table  8,  and  also 
their  variation  with  inclination  as  given  by  GEM  IOC.  Since  the  multiplying  factors  in 
equation  (14)  remain  large  up  to  a  higher  degree  than  those  in  equation  (12),  terms  up  to 
degree  53  in  GEM  IOC  are  included. 


The  variation  of  the  li 
GEM  10B/C 


?ed  harmonics  with  inclination,  and  implications  for 


Table  8 

_  _  I  1  _ 13 

Values  of  odd-degree  lumped  harmonics  <C , s ) 3 1  and  (CjS),,*  from  three  satellites 


Satellite 

I973-27B 

1968-103A 

1961  al 

i  (deg) 

50.0 

51.5 

97.2 

9-1.1 

10  C3, 

20800  +  10100 

-49000  ±  14000 

-47.3  ±  7.2 

q-1.1 

,0  S3. 

24700  +  11400 

11000  ±  7000 

-8.7  ±  4.2 

**31,31,15 

0.0002811 

0.0005269 

0.3430 

9-  -I.' 

10  F3!,31,15C31 

5.8  ±  2.8 

-25.8  ±  7.4 

-16.2  ±  2.5 

9_  -1.1 

10  F31 ,31 , 15S3l 

6.9  ±  3.2 

5.8  ±  3.7 

-3.0  ±  1.4 

9_-1.3 

,°  c3. 

36100  ±  4300 

10000  ±  3000 

63.9  ±  15.3 

9--1.3 

10  S3, 

17800  ±  2800 

2000  ±  3000 

42.0  ±  15.2 

F3 1 ,31 ,14 

0.001139 

0.001995 

0.2349 

9-  .-1.3 

10  F31 ,31 , 14C31 

41.1  ±  4.9 

19.95  +  6.0 

15.0  +  3.6 

9-  -"1.3 

in  f  s 

20,3  ±  3.2 

9.9  ±  3.6 

10  31,31,14  31 

*♦  »U  —  U  1 U 

It  is  not  possible  to  make  an  immediate  judgment  on  the  accuracy  of  GEM  IOC  from 
Figs  4  and  5:  an  indication  of  the  range  of  error  on  the  GEM  curves  is  required.  In  Tables 
9  and  10,  therefore,  the  values  of  the  lumped  coefficients  from  the  two  best  satellites, 
Samos  2  and  1973-27B,  are  compared  with  the  values  given  by  GEM  10C  at  the  appropriate 
inclination  (97.24°  or  50.04°),  on  the  assumption  that  the  individual  coefficients  in 
GEM  10C  have  sd  of  5  x  10 

Table  9 

Comparison  of  lumped  harmonics  from  Samos  2  and  GEM  IOC 


Lumped 

coefficient 


q-0.2 

10  C31 


gj.l 

10  C3. 


9--*. 3 

10  C3. 


GEM  10C 
for  i  ■  97 .24c 


-4.5  ± 
(-4.4  ±  6) 


-18  ±  8 
(-18  ±  8) 


-25  ±  12 
(-23  ±  12) 


Samos  2 


-2.9  ±1.2 


-47  ±  7 


64+15 


Lumped 

coefficient 


9-0.2 

1°  S3, 


g_— 1.3 

10S3I 


GEM  10C 
for  i  *  97 .24c 


-1.7  ±  6 
(-1.7  ±  6) 


-4  ±  8 
(-3  ±  8) 


-11  ± 
(-11  ± 


Samos  2 


9.0  ±  2.2 


42  ±  15 
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The  numbers  in  brackets  in  Table  9  are  the  values  that  would  be  obtained  by  using 
coefficients  of  degree  up  to  that  recommended  in  Ref  8,  namely  44,  55  and  61,  rather  .than 
the  maximum  degree  used  here,  namely  46,  53  and  53  respectively.  The  comparison  with 
GEM  IOC  made  in  Ref  8  utilised  a  preliminary  set  of  GEM  IOC  coefficients,  now  superseded: 
so  the  values  in  brackets  provide  a  revised  version  of  the  table  in  Ref  8  -  and  also  show 
that  altering  the  maximum  degree  does  not  have  much  effect. 


In  Table  9  the  GEM  IOC  values  of  C, 


-0.2 


and  S„ 


-0.2 


differ  from  the  rather  accurate 


■*31  31 

observational  values  obtained  from  Samos  2  by  0.3  and  1.8  times  the  GEM  sd.  This  suggests 
that  the  chosen  sd  for  GEM  is  nearly  correct.  For  the  odd  (C,S)  coefficients,  however, 
the  quoted  standard  deviations  from  the  two  sources  are  nearly  equal  and  their  sum  is  a 

,1,1  ,1,1  -“1,3 


'31 


and 


better  yardstick:  the  differences  between  ihe  GEM  10C  values  of  C,,  ,  S  ,  C 

__ 1 , 3  j i  j i 

S^|  and  the  corresponding  Samos  2  valuer  are  1.9,  0.4,  3.3  and  2.0  *  (the  sum  of  the 

sd) .  This  suggests  a  larger  sd  for  GEM  IOC.  In  the  odd-degree  lumped  harmonics,  the 

higher-degree  terms  are  much  more  important:  so  the  indication  is  that,  although  GEM  10B 

-9 

(degree  <36)  may  be  accurate  to  5  *  10  ,  GEM  !0C  is  probably  not.  For  Samos  2,  GEM  10B 

9  0,2  9-0, 2 

gives  10  C3|  *  -3.2  and  10  S3  -  -0.3  ,  and  these  are  both  closer  than  GEM  10C  to 


the  observed  values. 


Table  10 

Comparison  of  lumped  harmonica  from  1973-27B  and  GEM  IOC 


Lumped 

coefficient 

GEM  10C 
for  i  ■ 50.04° 

1973-27B 

Lumped 

coefficient 

GEM  IOC 
for  i  - 50.04° 

197J-27B 

,  0,2 

10  C3. 

-0.7  ±  0.7 

0.86  ±  0.32 

,  0,2 

1°  s3. 

0.1  ±  0.7 

-1.21  +  0.20 

10  Si 

OO 

+1 

r* 

1 

21  ±  10 

gj.l 

1°  S3. 

-8+8 

25+11 

10  Si 

-1  +  2 

-36  ±  4 

,  -1,3 
10°S 

31 

-2  +  2 

18+3 

_0, 2  _0,2 

A  similar  comparison  for  Table  10  shows  that,  for  C^  and  S3)  ,  the  GEM  10C 

values  differ  from  I973-27B  by  2.2  and  1.9  standard  deviations.  (The  inclusion  of  terms 

6  0,2 

up  to  degree  60  in  GEM  10C  gives  no  significant  improvement:  10DC,!  changes  to  -0.4 
5_0,2  _  _  1  1 

and  10  S.,  to  0.6.)  For  (C,S)  J  the  differences  are  1.6  and  1.7  »  (the  sum  of  the 
Jl  -  -  -I  3  _  J 1 

sd) ;  and  for  (C.S)^’  the  differences  are  5.8  and  4.0  x  (the  sum  of  the  sd) .  This  can 
only  be  called  serious  disagreement:  however,  since  all  the  lumped  harmonics  involve  high- 
degree  coefficients,  the  previous  conclusion  about  GEM  10B  (degree  <36)  remains  unaffected. 

Thus,  although  there  remains  room  for  doubt,  our  provisional  conclusion  for  31st- 

order  coefficients  is  that  GEM  10C  cannot  be  recomnended ,  but  the  GEM  10B  coefficients 

-q 

(degree  <36)  are  probably  within  the  5  *  10  accuracy. 


i. 
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6  SUMMARY  OF  RESULTS 

6. 1  Harmonica  of  order  30 

The  values  of  even-degree  30th-order  lumped  harmonics  obtained  from  analyses  of 
ISth-order  resonance  are  numerous  enough  to  allow  a  solution  for  individunl  coefficients 
(Table  3).  The  best-determined  coefficients  are  those  of  degree  30,  namely: 

10,S0,30  ‘  '  '-2  *  l0%0.30  •  9'6±  U3  ' 

These  standard  deviations  correspond  to  an  accuracy  of  about  I  cm  in  geoid  height. 

The  variation  of  the  even-degree  lumped  harmonics  with  inclination  given  by  the 

Goddard  Garth  Model  10B  is  generally  similar  to  the  variation  given  by  the  resonance 

-9 

solution  (see  Fig  1),  and  it  appears  that  a  standard  deviation  as  low  as  3  *  10  may  be 

appropriate  for  the  30th-order  coefficients  in  GEM  I0B.  For  the  (30,30)  coefficients, 

-9  -9 

the  values  from  GEM  I 0B  (with  sd  3  *  10  ),  and  from  Rapp's  1977  solution  (sd  5  *  10  ), 

are  consistent  with  those  quoted  above. 

At  the  single  inclination  where  a  test  is  possible,  the  odd-degree  30th-order 
S-coeff icients  from  GEM  agree  well  but  the  C-coef ficients  seem  to  be  too  small  (see  Fig  2). 

All  in  all,  progress  in  evaluating  30th-order  coefficients  is  encouraging,  though 
more  analyses  of  long-enduring  15th-order  resonances  are  required. 

6.2  Harmonics  of  orders  29  anJ  31 

The  picture  is  different  for  these  harmonic.?.  There  are  not  enough  analyses  of 
29:2  or  31:2  resonances  to  allow  the  determination  of  individual  coefficients,  and  all 
that  can  be  done  is  to  compute  the  variation  of  lumped  harmonics  wich  inclination  as  given 
by  GEM  10C  (GEM  I0B  only  goes  to  degree  36,  and  higher  coefficients  are  usually  needed), 
and  to  compare  lumped  harmonics  from  GEM  10C  with  the  small  number  of  accurate  values 
found  from  resonance  analyses.  Figs  3  to  5  show  that  there  is  no  consistent  agreement. 

For  3l8t  order,  there  is  an  indication  that  the  GEM  1 0B  values  (degree  31  to  36)  may  be 
accurate  to  5  *  10  ,  but  the  GEM  10C  values  (degree  >36)  disagree.  Since  GEM  10C  was 

derived  as  a  model  of  the  ocean  surface  to  degree  and  order  180,  the  use  of  a  small  group 
of  its  coefficients,  of  degree  37  to  50,  say,  may  be  inappropriate.  For  29th  order,  the 
resonance  results  are  fewer  and  less  accurate,  so  the  disagreement  is  inconclusive. 

Further  analyses  of  29:2  and  31:2  resonance,  particularly  the  former,  are  needed  if 
future  comprehensive  models  of  the  grc-.^ty  field  are  to  be  tested  by  the  resonance  method. 

6.3  Conclusions 

Resonance-testing  of  comprehensive  gravitational  field  models  has  previously  been 
applied  (ej  Refs  3  and  26)  for  various  orders  up  to  15.  Here  the  test  is  made  for  orders 
29,  30  and  31,  with  mixed  results:  GEM  10B/C  emerges  well  for  order  30;  for  order  31,  the 
GEM  10B  values  (degree  .£36)  are  probably  good,  but  the  GEM  IOC  values  (degree  >36)  are 
probably  not;  for  order  29  the  test  is  indecisive. 
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109Cr 


-  1 .2  ±  1.1 


9- 

10*S, 


-  9.6  ±  1.3 


30,30  30,30 

The  standard  deviations  here  are  equivalent  to  1  cm  in  geoid  height. 

For  the  29th-  and  31st-order  harmonics,  and  for  the  30th-order  harmonics  of  odd 
degree,  there  are  not  enough  values  to  determine  individual  coefficients,  but  the 
lumped  values  from  particular  satellites  can  be  used  for  'resonance  testing'  of 
gravity  field  models,  particularly  the  Goddard  Earth  Model  10B. 
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